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Abstract 


Some problems on variations are raised for classical discrete mechanics and 
field theory and the difference variational approach with variable step-length 
is proposed motivated by Lee’s approach to discrete mechanics and the differ¬ 
ence discrete variational principle for difference discrete mechanics and field 
theory on regular lattice. Based upon Hamilton’s principle for the vertical 
variations and double operation of vertical exterior differential on action, it is 
shown that for both continuous and variable step-length difference cases there 
exists the nontrivial Euler-Lagrange cohomology as well as the necessary and 
sufficient condition for symplectic/multi-symplectic structure preserving prop¬ 
erties is the relevant Euler-Lagrange 1-form is closed in both continuous and 
difference classical mechanics and field theory. While the horizontal variations 
give rise to the relevant identities or relations of the Euler-Lagrange equation 
and conservation law of the energy/energy-momentum tensor for continuous 
or discrete systems. The total variations are also discussed. Especially, for 
those discrete cases the variable step-length of the difference is determined by 
the relation between the Euler-Lagrange equation and conservation law of the 
energy/energy-momentum tensor. In addition, this approach together with 
difference version of the Euler-Lagrange cohomology can be applied not only 
to discrete Lagrangian formalism but also to the Hamiltonian formalism for 
difference mechanics and field theory. 
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I. INTRODUCTION 


Variation problems play a fundamental, even central, role in various types of continu¬ 
ous mechanics (see, for example, P], 0, |]l|) and field theories for both Lagrangian and 
Hamiltonian formalisms, which can be transformed each other in many cases via Legendre 
transformation. A lot of important issues such as the equations of motions, the (intrinsic) 
symplectic or multisymplectic preserving properties, conservation laws associated with cer¬ 
tain symmetries, topological properties etc. are very closely related to the variation problems 
with hxed or variable domains. On the other hand, however, there are some open questions 
relevant to variation problems in different kinds of discrete mechanics, held theories as well as 
corresponding symplectic and multisymplectic algorithms, although discrete variation prob¬ 
lems still play an important role, particularly, for the discrete Lagrangian formalism of these 
discrete systems. In fact, for a long period, there had been no any discrete variation ap¬ 
proach available to both discrete Lagrangian and Hamiltonian formalisms that relate each 
other by discrete Legendre transformation until the difference discrete variational approach 
has been proposed very recently |^, [|H . 


As far as discretized version with difference for mechanics is concerned, Lee suggested 
a discrete variational approach to discrete Lagrangian mechanics and relevant algorithm in 
early 1980’s ii, ||20|| , ||21|| . In Lee’s approach the time steps are variable in view of time 


being a dynamical variable in order to preserve the energy discretely. Veselov 03, 03 also 


proposed the discrete variational principle by the end of 1980’s, which is almost the same 
as Lee’s approach except without taking variation with respect to the discrete time so that 
it does not keep the conservation of the energy discretely in general. In addition, both 
approaches are merely available to discrete Lagrangian mechanics and nothing to do with 
discrete Hamiltonian mechanics. 

On the other hand, Ruth and Feng proposed the symplectic algorithm for Hamil¬ 
tonian mechanics. In this algorithm (for a review, see [^), the time step-length is hxed 
and the symplectic preserving property is discretely kept. However, the discrete version of 
energy conservation can not be maintained discretely in general. The symplectic algorithm 
plays very important role in computational mathematics and its applications cover various 
branches in sciences. It is well known that the progresses of symplectic algorithm promote 
further development of the structure-preserving algorithms. In these algorithms, the time 
step-length is always being hxed so that the price paid by keeping structure preserving is 
the loss of other conservation laws for the continuous cases in general. 

Structure preserving criterion: It is widely accepted that the discrete systems should 
be thought as the discrete counterparts of the corresponding continuous systems. However, 
in order to discretize continuous systems, a guide line in needed. In Feng’s hrst paper on 
the symplectic algorithm, he wrote a working hypothesis: “It is natural to look forward 
to those discrete systems which preserve as much as possible the intrinsic properties of 
the continuous system.” [|T3 In fact, this statement should be regarded as a criterion, the 


structure-preserving criterion, for constructing mostly quarried one in all kinds of structure¬ 
preserving algorithms. However, in order to carry through this criterion, it is needed to know 
how to answer the following simple questions. 

Problem 1: What are the most important intrinsic “structures” in continuous systems, 
such as classical mechanics and held theory, that should also be maintained in the course of 
discretization? What are the discrete counterparts of these “structures” and how to preserve 
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them in certain discrete version? What about the lowest price has to be paid in the course 
of discretization? 

It is well known that there are two classes of conservation laws in canonical conservative 
mechanics. The hrst class of conservation laws is that of phase-area conservation laws char¬ 
acterized by the symplectic preserving property and another class is related to energy and 
all hrst integrals of the canonical equations. Thus, the following questions can apparently 
be raised. 

Problem 2; Is it possible to establish a kind of discrete mechanics and/or structure¬ 
preserving algorithms in such a way that they not only discretely preserve the symplectic 
property but also some other conservation laws, specially the energy conservation? Do 
these discrete systems can be established by a discrete variational approach? Does this 
discrete variational approach can be applied to both discrete Lagrangian and Hamiltonian 
formalisms? 

Although in Lee’s approach, it should be able to prove that in addition to discrete energy 
conservation the symplectic structure is also preserved since Lee’s approach is a discrete 
variational approach But the framework of either Lee’s approach or Veselov’s one can not 
be applied to the discrete Hamiltonian systems. To our knowledge, therefore, these problems 
are still open, at least partly. 

In the cases of discrete held theories and multisymplectic algorithms, the symplectic 
algorithm has been generalized to the multisymplectic one in what is called “Hamiltonian 
formalism” P|. On the other hand, Veselov’s discrete variation approach to the discrete 
mechanics has also been generalized to held theories in Lagrangian formalism to get so-called 
“ variational multisymplectic integrators” |^, [Q. In both approaches, the step-lengths 


are hxed so that the energy-momentum tensor cannot be conserved in general although the 
multisymplectic structure preserving property in held theory can be maintained discretely 
in certain manner. Thus a set of similar questions can also be raised to the discrete held 
theory and multisymplectic algorithms. 

Problem 3; Is it possible to establish a discrete variational approach to describe a 
kind of discrete held theories and/or multisymplectic algorithms in such a way that not 
only the multisymplectic property is discretely preserved but also the conservation laws such 
as energy-momentum conservation law are discretely maintained in certain version? Is it 
possible to apply such an approach in both discrete Lagrangian and Hamiltonian formalisms? 
To our knowledge, these problems are still open as well. 

As was just mentioned, for the diherence discrete variational approach to discrete mechan¬ 
ics and held theory with hxed step-length diherences, it can be applied to both Lagrangian 

0 , 


I^. The key point of this approach is to regard the diher- 


and Hamiltonian formalisms 
ences with hxed step-length as a kind of entire geometric objects, which play an analogical 
role with the one played by derivative in continuous cases. 

Problem 4: Is it possible to generalize the diherence discrete variational approach with 
hxed step-length diherences to the one with varied step-length diherences so as to the discrete 


^In 1^] this problem has been partly solved by defining a conserved discrete energy in Veselov’s 
approach. The complete resolving to this problem in Lee’s framework has been made until very 


recently in i, 0- 
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energy conservation law in certain version may be kept together with the symplectic and/or 
multisymplectic preserving properties? 

In order to carry through the structure-preserving criterion, we present an approach to 
the discrete variation problems. Namely, the difference variational approach with variable 
step-lengths named variable difference variational approach to these problems in discrete 
mechanics and field theory. In other wards, this approach is available to the discrete to¬ 
tal variations with keeping the step-length of differences to be varied by an equation from 
the variation of the discrete action with respect to discrete time and/or discrete space co¬ 
ordinates. As was mentioned, this approach is a natural generalization of the difference 
variational approach with fixed step-lengths proposed recently in |^, |]T^ for the ordinary 
discrete variation problems with fixed discrete domain. In fact, the approach in [|^, |]^ is 
just discrete vertical variation so that it is natural to keep the step-length being fixed. The 
most important key point of this approach is that in the course of calculation of variation 
problems in discrete mechanics and field theory, the differences with variable step-lengths 
are kept as entire geometric objects as much as possible. Consequently, this discrete varia¬ 
tion approach not only keep the advantage of Lee’s discrete variation, which conserves the 
energy of the system discretely, but also the advantage in variation in Veselov type, which is 
symplectic or multisymplectic. In addition, this variable difference variation approach can 
be applied not only to the Lagrangian formalism but the Hamiltonian formalism as well for 
both discrete mechanics and discrete field theory, since the discrete canonical “ momenta” 
and discrete version of Legendre transformation can be introduced in terms of variable step- 
length differences. For simplicity, we consider in this paper the discrete Lagrangian of first 
order of difference only. 

The paper is organized as follows. In section 2, we recall the total variation problems 
in Lagrangian and Hamiltonian mechanics. In section 3, we present the variable difference 
variational approach and deal with the difference Lagrangian and Hamiltonian mechanics. 
In section 4, we recall briefly the total variation problems in Lagrangian and Hamiltonian 
field theory with generic variables. In section 5, we apply the variable difference discrete 
variational approach to the total discrete variation problems in difference discrete Lagrangian 
and Hamiltonian field theory with generic variables. Finally, we end with some remarks. 


II. GENERAL VARIATIONS FOR CLASSICAL MECHANICS 

Let us recall briefly the general or total variation calculus with variable domain in classical 
mechanics. 

Let t G T ~ i? be the time, M an n-dimensional configuration space. Consider a fibre 
bundle E{T, Q, vr) with projection tt : i? —>■ T on T, 7r“^ : t ^ Qt isomorphic to M is the 
fibre on f G T. Denote F(i?) the sections on E, TE the tangent bundle of E, TyE cTE the 
vertical sub-bundle of TE, etc.. 


A. Variations in Lagrangian mechanics 

We first consider the Lagrangian mechanics. The Lagrangian of the system is denoted 
as L{q‘^{t),q^{t)]t), i = which is a mapping from TE to R. For simplicity, the 

Lagrangian is of the first order. The action functional is 
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( 2 . 1 ) 


Jtl 


Here g*’s are coordinates in the fiber, g*(t) describes a curve with ending points a and 
b, ta = = t 2 , along which the motion of the system is assumed to take place, and 

= dq''{t)/dt. 

Let us consider the general variation of q'‘{t) 

q\t) ^ q'\t') = q\t) + 5tq\t) (2.2) 

accompanied with an infinitesimal re-parameterization of time t 

t —>• t'{t) = t + 5t. (2.3) 

Here Stq\t) denotes the total variation that can be divided into two parts, 

btq'{t) = 6yq^{t) + 6hq'{t), 
bvqKt) = q'Kt)-qKt), 
dkqKt) = q'\t')-q'Kt) 

= q\t') - q\t) + 0(5^) = 5tiq\t) + 0(52), 

where 6yq^{t) denotes the equal time part variation or the vertical one and 6hq^{t) the hor¬ 
izontal part along the fibre induced by the re-parameterization of the time t (|2.3|) . If we 
introduce a variational vector field on T 


5(i):=*L (2,5) 

the horizontal variation 6hq'‘{t) is the Lie derivative of g*(f) with respect to the variation 
vector field (p75|). 

Similarly, 

btq'it) = {Sv + Sh)q''{t), 
bvq\t) = q'\t)-q\t), 
bhq\t) = i,q'\t') - ^q'\t) 

= + 0 ( 5 ^) = + 0 ( 5 ^). 


Note that 6hq^{t) is also the Lie derivative of g*(f) with respect to the variational vector field 
(P^.SI). In fact, this is true for a kind of functionals of q^{t), q^{t) and t: 


ShF(q'(t),q‘(t), t) = L(F(q'(t), q'(t),t). 


(2.7) 


To the time change 
formula 


is associated the change in the measure in (|2.1|) given by the Jacobi 


, / db . /. d ^ , 

dt = —dt = (1 -|- —5t)dt, 
ot dt 

i.e. 

6{dt) = d{t + 6t) — dt = dt—6t. 


( 2 . 8 ) 


(2.9) 
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It is easy to see that this change in the measure is also the Lie derivative of the measure 
with respect to the variation vector (PTB|): 

5{dt) = L^dt = dL^t = d5t. 


Here the commuting property between d on T and the Lie derivative is used. 
Now the Lagrangian is changed as 

LW(t ),!«'((); t) - r(9'‘(*'). *') = t) + s,l, 


and the action is also deformed to 


= Ju %dt{L{q\t), + 6tL} 

= jf^dt{L + {i5t)L + 5,L} 

= S{[q\t)]-,ti,t2) +6tS 


A more or less straightforward calculation shows 

r ^2 rl f) d f)T 

S,S = I + IjH + -m + -(^V - HM)}, 

where [Lgi] is the Euler-Lagrange operator and H is the energy (Hamiltonian) 

rr 1 _ d^dL Tj -i T 

^ Siji dt^dq‘^' dq''^ 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


The vertical and horizontal variations should be separated as the independent ones and 
this leads to 


pl 2 rj f) r 

S,S = dt{lL,.]6„q‘ + -(—S,qi)}, 


(2.14) 


and 

pt2 ^ f) r! FiT 

= I dt{lL,]S.q' + Ij^H + -m + -(^4,^- - mt)}. (2.15) 

For the vertical part, the Hamilton’s principle leads to the Euler-Lagrange equation if 
dvQ^lti dyq^\t2 0 


dq^ 


d dL. 
dt dq^' 


= 0 . 


For the horizontal part, however, it is easy to check 

\LM + iH + §;L = 0, 

■d := §^Shq^ — H6t = L6t. 


(2.16) 


(2.17) 


Therefore, 
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(2.18) 


4S = dt^(LSt) = 0. 

Jti at 


This is just the invariance of the action S with respect to the re-parameterization of time. Of 
course, from the hrst equation of (|2.17|) , it still follows the conservation law for the energy if 
and only if the Euler-Lagrange equation is satished and L does not depend on t manifestly. 
If the identities in (|2.17|) is employed directly, it follows (see, for example, 0) that 


pf2 r} f)T 

6^S = dt{[L,.]6,q^ + + L6t)}. 


A quantity now can be dehned 


^ dL ■ 

J := -^Syq^ + L6t, 
oq^ 


(2.19) 


( 2 . 20 ) 


the invariance of S under the re-parameterization of t, i.e. 6tS = 0, leads to the conservation 
of the quantity J if and only if the Euler-Lagrange equation is satished: 


d 

dt 


J = 0. 


( 2 . 21 ) 


Remark 2.1: 

Introducing an exterior differential operator dy along the hbre that satishes 

dy^ = 0, {dy,dh} = 0, d:=dy + dh, ( 2 . 22 ) 

where dh and d is the nilpotent exterior differential operator on T*T and T*E respectively 
as was as a vertical variational vector held 

:= '5.9‘(i)^. (2.23) 

then 


Syq\t) = i^^dyq" = i^^dq\ 


(2.24) 


By means of the vertical variational vector held (|2.23|) on TQ, 5yS can also be expressed 
as its contraction with 1-form dyS G T*Q 


kgdyS = SyS. 


(2.25) 


We may calculate dyS G T*Q. Since dy commutes with the integral of dt (see also, for 
example, the functional diherential calculus in |^), it is straightforward to get 


dyS = J'" dt{[Lg.]dyq^ + 


where 6 is the Lagrange 1-form 


6 ^.^dyq\ 


(2.26) 


(2.27) 
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Now by contracting with the vertical variational vector held (|2.23|) it follows straightforwardly 
6,S in ( CT )- 

Furthermore, due to the nilpotency of dy, it is easy to get 


, ^ d 

Cly£ — 0, 

at 


where S is called the Euler-Lagrange 1-form [jT^, |T^, |]T^, dehned by 

:= [Lqi]dyq\ 

10 is the symplectic 2-form and in local coordinates: 

d^L d^L 

u; := dyO = -K—^dyq^ A -F ^.-^.- dyq^ A dyq\ 
oqWq^ oqWcf 


(2.28) 


(2.29) 


(2.30) 


From (|2.29|) , (|2.26|) and (p.28|) , the following theorem can be proved |T^, [^], [|T4 


Theorem 1: For all Lagrangian of a kind of systems with first order of derivatives on the 
bundle E{T, Q, tt) 

1. The following Euler-Lagrange cohomology is nontrivial: 


Hlm ■= {£\dy£ = t]}/{£\£ = dya), 


where a = a(g*(t), g*(t); t) is an arbitrary function of {q\t),q\t);t). 

2. The necessary and sufficient condition for conservation of the symplectic 2-form, i. e. 

j^io = 0, (2.31) 

is the corresponding Euler-Lagrange 1-form being closed. 

Remark 2.2: 

From the dehnition of the Lie derivative it can be seen that the horizontal variations are 
given by the Lie derivative with respect to the variational vector held. 

Let be a vector held on T, exp{Xf) be the how with parameter A, i.e. the one-parameter 
diheomorphism group, induced by ^, / a diherential or a vector on T. The inhnitesimal 
change of / under how is described by its Lie derivative with respect to the vector held f 


Lifif) ■= lim ^{01/(exp(AOt) - f{t)} = ^|A=o(0I/(t')), t' = exp{\f)t. (2.32) 

Here (f)\ is the bull-back or the inverse diherential for the diherential form or vector respec¬ 
tively. 

Taking i = Hf) = in (|2.5|), it follows that the Lie derivative of fit) with respect to ft 
gives rise to the horizontal variation of f{t). 

On the other hand, the time variation 5t can be expressed by the contraction between 
the variational vector held (|2.5| ) and 1-form dht on T*T, i.e. dt on T*E 


i^^dht — 6t. 


( 2 . 33 ) 









































It is also feasible to express the variation 6hq^{t) as contraction of a horizontal variation 
vector held with or dq\ To this purpose, ^h,q along the hbre with respect to horizontal 
variations of 


^h,q ■= Shq'{t) 


A 

dq^ 


(2.34) 


should be introduced. Combining with the vector held in (|2.5|) , the general horizontal 
variational vector held should be dehned as 


'■= 6 + ^h,q = (2.35) 

Its contraction with d^q^ or dg* leads to 

i^^dq" = dq" • = dhq\t). (2.36) 

In general, for any functional of g*(f) and g*(f), F{q\t),q\t)) : TQ —> R, its (horizontal) 
variation induced by ( p.3D is 


F{q^{t),q%t)) ^ F{q\t'),^q^{t')) = F(g*(f), g*(f)) + 4F(g*(t), g*(f)), 
5hF{q\t),q\t)) = i^JF{q\t),q\t)). 


(2.37) 


Remark 2.3: 

For the total variation, a total variational vector held for g*(f) along the hbre can also 
be introduced 

itotal '■= ivF ih = 

= + Mt) + ( 2 . 38 ) 

whose contraction with dq^ leads to the total variation dtq^{t) 

ktotaM = dq" ■ itotai = dtq\t). (2.39) 

If we introduce the Lagrangian 1-from 

L := L(g*, g\ f)(if (2.40) 

and take 0 = d^L, it is easy to see that the theorem 1 still holds. This means that the 
total variations keep the Euler-Lagrange cohomology as well as the necessary and sufficient 
condition for symplectic structure preserving property in classical mechanics. 

Remark 2.4: 

In some literatures (see, for example, @), it is required that Hamilton’s principle holds 
for the total variation of the action, i.e. 5tS = 0, and regard djg* and 5t as independent 
variations. Thus it follows the Euler-Lagrange equation, the conservation relation for the 
energy and the surface term 
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(2.41) 


dL 


d { dL' 


dq^ dt f ^ ’ 

^ - EH. 


If -^L = 0, i.e. the system is conservative, the energy H is conserved. However, 5tq^ is 
actually dependent on 6t. Therefore, it would be better to regard and 6t as independent 
variations. 


B. Variations in Hamiltonian mechanics 


The action principle should, of course, be applied to the Hamiltonian mechanics. In order 
to transfer to the Hamiltonian formalism, we introduce a set of conjugate momenta from the 
Lagrangian L{q\t),q\t)-,t) 

dL 

~ w 

and take a Legendre transformation to get the Hamiltonian 

H := H{q\pj;t) = Pkq^ - L{q\ q^;t). 

Now the action functional can be expressed as 


(2.42) 

(2.43) 

(2.44) 


S{[Pi{t)],[q"{t)];ti,t 2 ) = f dt{pkq’" - H{q\pj-,t)}. 

J t\ 

The total variation of the action can be calculated 
5tS = 5^S + 5hS, 

^vS = dt{[Hp^]5^Pi - [Hqi\5^q^ + ji^pAq")}, 

hS = Hi dt{[Hpl^5hP^ - [HHdhq^ + [^H - + Hp^S^q^ - H6t)}, 

where [iLpJ, [Hqi] are canonical operators 

f)J-T f)M 

Thus, the stationary requirement for the vertical variation of the action S^S = 0 leads to 
the canonical equations 

dH . dH 

^ ■ 

While the time re-parameterization invariance of the action, i.e. 6hS = 0 gives rise to an 
identity on the condition of the energy 


(2.45) 


(2.46) 




(2.48) 


and the boundary term that leads to so called “extended symplectic potential” is in fact a 
total divergence 

r dt^Mnq^ - HSt) = ^ dt^iLdt) = 0. 

Jti at Jti at 

Similar to the Lagrangian mechanics, all remarks in last subsection can be made for the 
Hamiltonian formalism. Especially, the theorem 1 can also be established here. 


(2.49) 
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III. GENERAL VARIATIONS FOR DISCRETE MECHANICS 


We have proposed a difference variational principle for the (vertical) variation in discrete 
Lagrangian mechanics in view of the differences with fixed time step-length being regarded 
as entire variables [Q, |]^. One of advantages of this approach is that it is available to 
both Lagrangian and Hamiltonian formalism for the discrete mechanics. This approach can 
also be generalized to the total variation for the differences with variable time step-length 
being regarded as entire variables in such a way that the variable time step-length should 
be determined by an equation given by the variation problem with variable discrete integral 
domain. In Lee’s wards, discrete time is regarded as a dynamical variable. 

Consider the case that “time” t is difference discretized 


t E R ^ t E Tf) — {{tk, tfc+i — tk k E Z)'\ (3.1) 

and the step-lengths are determined by a variational equation, while the n-dimensional 
configuration space Mk at each moment tk, k E Z , is still continuous and smooth enough. 

Let N be the set of all nodes on To with index set Ind{N) = Z, M = Ufcez Mk the 
configuration space on Td that is at least pierce wisely smooth enough. At the moment tk, 
Mk be the set of nodes neighboring to tk- Let R the index set of nodes of A4 including tk- 
The coordinates of Mk are denoted by q''{tk) = = 1, • • •, n. T{Mk) the tangent bundle 

of Mk in the sense that difference at tk is its base, T*{Mk) its dual. Let Aik = Uze 4 Afz be the 
union of configuration spaces Mi at ti, I E R on Mk, TAik = Uzg 4 TMi the union of tangent 
bundles on Mk, F{TMk) and F{TAtk) fhe function spaces on each of them respectively, 
etc.. In the difference variational approach, we will use these notions. 


A. Variable difference Lagrangian mechanics 


Let us consider the system with a discrete Lagrangian on F{T{Mk x Td))- For 

simplicity, the Lagrangian is of the first order of differences 


with the difference of at tk defined by 

;= 1 - 1 -. 

4+1 - 4 

The discrete action of the system is given by 


(3,2) 


(3.3) 


Sd = Y.itk+i-tk)LD^^\q^^^\— - ^^;tk). (3.4) 

kez tk 

The discrete total variations for = g*(fzc) should be defined as follows 

5zg'(") := q’Mk) - qKtk) = 4 = 4 + 54, 

;= g'*(4) - g*(4), 5,;4 = 0, (3.5) 

5^gi(Zc) _ ^ ^ 


11 














(3.6) 


It can be shown that horizontal variation is given by 

The discrete total variations for are dehned as 




g'h*fc+i)-g'h*fc) 


^fc+l 


^fc+l 


(3.7) 


Dne to the dehnition of the difference with variable time step-length (|3.3|) and the Leibniz 
law for it 


= (A,/(^))gW + E/W(A,gW), 

where E is the shift operator dehned as 

j{k) _ j{k+l) ^ £-1 j{k) _ ^{k-1) ^ 


it follows 


Namely, 


6tAkq^^^^ = Afc((54g*(^)) - {Ak6tk)Akq^^^\ 
4Afcg*(^) = 54+iA(Afcg*(^)). 


6,Akq^^'^'> = Ak5,q^^^\ 

^Afcghfc) = Afc^/^ghfc) _ (Afc(54)Afcghfc). 


Using above properties and 

4) ^A:(<5i4) (4+1 ^k)-! 

the total variations of the discrete Lagrangian can be calculated as follows 


-h Ak{pi^^^Akq^^^-^'^)5tk + 

+ Afc(p,("+i)<5tg*W -p,WAfcg*("-i)54), 

where [L^^k)] is the discrete Euler-Lagrange operator 

• Q^i(k) '^aAghfe-i)^’ 
and the discrete canonical conjugate momenta 

p,T (fc-i) 

„.(U .= 

* dAq'^g^~^'> 

Thus the total variation of action is given by 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 
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(3.16) 


StSn = Ek{tk+i-h){{ASt,)Lo^’^^ + StLn^^'>} 

= Efc(4+i - h){[L^.i.)]6tq^^^^ + + ^)54 

where is the difference Hamiltonian that can be introduced through the discrete Leg¬ 

endre transformation 




(3.17) 


Thus the total variation of the discrete action (p.l6|) can be written as 


^tSo — ^vSd + ^hSoi 

SvSd = Efc(4+i - 4)[Lgi(fc)](5^g*(^)-h A(pi(*^+L5^g*W), 

ShSd = Ekitk+i - h){[L^m]hq^^’^^ + + ^)54 


(3.18) 


The variational principle requires SySo = 0 and the discretized re-parameterization in¬ 
variance with respect to discrete time may also leads to ShSr, = 0 if this invariance does 
exist. Thus it follows the discrete Euler-Lagrange equations for gh^)’g 


dl^ 

dq^{k) 


A( 




= 0 , 


and the equation for the variable time step-length 


(3.19) 


( 


dqi-{k) 


A( 




))Ag*W + AHd^’'~^^ 


dH^ 

dtk 


0 . 


(3.20) 


It is more or less straightforward to show that if the time step-length is fixed the equation 
( p.20| ) has no solution in general even if the Lagrangian does not depend on discrete time 
manifestly. In other wards, for the conservative discrete Lagrangian mechanics the time step- 
length should be variable in general so as to the energy of the system can be kept conserved 
discretely. 


Remark 3.1: 

We may introduce exterior differential operators d, dy and du on T*{M x Tp), T*M and 
T*Tp respectively. They are nilpotent and satisfy 


d = dy + dh, {dy,dh} = 0. (3.21) 

Especially, dh is due to the difference on Tp and satisfy Leibniz’s law for ordinary forms 0. 

Remark 3.2: 


^It is needed some noncommutative differential calcnlns to completely clarify the properties of dh- 
For the case that At is fixed, the noncommutative differential calculus can be found in [15|, [16|. 
For the case of variable time steps, similar noncommutative differential calculus can be established 

0 . 
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Actually, analog to the case with hxed time steps [Q, it can be established the dif¬ 


ference version for the Euler-Lagrange cohomology and the necessary and sufficient condition 
for the difference conservation law of the discrete symplectic 2-form. 

From 6^3D in (|3.18|), it is easy to see that we may take dy on Sd to get 


dyS £) ^ + 1 tk)dyL £) 

k 


(k) 


dyLn^’^^ = 


(3.22) 


where Sd^^\ are the discrete Euler-Lagrange 1-form and symplectic potential 1-form 

respectively 


:= [L^m]dyq^^’'\ := dyq^^'^\ (3.23) 

Then due to the nilpotency of dy, it is straightforward to get 

dySo^'’^^ + = 0, := dyOn’'’^^ = A dyq^^'^K (3.24) 


Therefore, we may get the discrete version for the theorem 1 [^, |T^, [Q: 


Theorem 2: For all diserete Lagrangian of a kind of systems with first order differences 
on the bundle EiT^i, Q,7i) ^ M x T^, 

1. The following discrete version of the Euler-Lagrange cohomology is nontrivial: 

Hdcm:={C losed Euler-Lagrange forms}/ {Exact Euler-Lagrange forms}. 

2. The necessary and sufficient condition for conservation of the discrete symplectic 
2-form, i.e. 


AkUJD^'^^ = 0, (3.25) 

is the corresponding discrete Euler-Lagrange 1-form being closed. 

Remark 3.3: 

In this paper, Td is an inhnite chain. It is reasonable to consider an interval on To. We 
will publish the issues on this topic elsewhere. 


B. Variable difference Hamiltonian mechanics 


Now we consider the total difference variation on the phase space in the discrete Hamil¬ 
tonian formalism with variable (time step-length) difference. 

In order to transfer to the discrete Hamiltonian formalism, it is needed to introduce 
the discrete canonical conjugate momenta according to the equation ( p.l5| ) and express the 
discrete Lagrangian by the discrete Hamiltonian via Legentre transformation (|3.17| ). Thus, 
the discrete action can be expressed as 


Sd 


Ek/k+i - t,)Ln^^\q^(^\ 
Ekitk+i - tk){p/^^+LA^/(k) _ 


(3.26) 


And its total variation reads 
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(3.27) 


StSo — SySn + ShSo 
= Efc(4+i - 

-(Ap,W + + (Ai7^(^-') + ^)54 

Similar to the discrete Lagrangian formalism, Hamilton’s principle requires 6^3d = 0 and 
the discretized re-parameterization invariance with respect to discrete time may also lead to 
ShSn = 0 if this invariance does exist. Thus it follows the discrete canonical equations for 
pg^^’s and 


Aghfe) = 


OH, 


D 


(fc) 




Apg'^) = 


dH^ 


D 


(k) 


dq^{k) 


dtk 


(3.28) 


0. (3.29) 


the equation for the variable time step-length 

- (Apg") + ^4^)Ag*(") + 

It is also more or less straightforward to show that if the time step-length is fixed the 
equation (|3.29|) has no solution in general even if the Hamiltonian does not depend on discrete 
time manifestly. In other wards, for the conservative discrete Hamiltonian mechanics the 
time step-length should be variable so as to the energy of the system can be kept conserved 
discretely. In [|^, 


11 , 18 , this issue has been studied. 


It should be mentioned that all remarks in last subsection can be made here and the 
theorem 2 can also be established for the discrete Hamiltonian formalism. 


IV. GENERAL VARIATIONS FOR FIELD THEORY 

Consider a bundle E{X, Q, vr), the fibre Q ~ M. For simplicity, let X = be an n- 

dimensional Minkowskian space as base manifold with coordinates x^, (p. = 0, • ■ •, n — 1), M 
the conhguration space on with a set of generic (scalar) helds u^{x), {a = 1, ■ • ■, r), 

TM the tangent bundle of M with coordinates (m",m“), where = |^, F{TM) the 
function space on TM etc. 

We also assume these helds to be free of constraints. In fact, the approach here can easily 
be applied to other cases. 


A. General variation in Lagrangian formalism 


The Lagrangian of the theory is supposed to be the hrst order of derivatives of the helds 
and dependent to the coordinates manifestly, i.e. L(m", x^), and the action is 

5([«“(a:)];a:^) = J^d^xL{u^,u^-,x^). (4.1) 

Let us consider the variations of the helds, i.e. total variation Stu°‘, vertical one and 
horizontal one ShU°‘'- 
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(4.2) 


m " —>• u'°'{x’) = u°‘{x) + 6 tu°‘{x), 

6tU°‘ = 5yU°' + 

5 -vU°‘{x) := u'°'{x) — u°'{x), 

5 hU°‘{x) := u'^\x') — u'°‘{x) = u°‘{x') — u°'{x) + 0 {S'^) 

= 6x^di^u°‘{x), 

accompanying with the coordinates’ inhnitesimal continuous transformation 


x^ —^ x'^ = x'^ + Sx^. 


The corresponding changes in the derivative of helds u" are 

6t{^u°{x)) = 6 ^{£j:U^{x)) + 4(^M“(a;)), 

^y Aw^ix)) := 

hi^u^ix)) = dx^£^{£jru^{x)). 

Now the action m is also changed as follows 


S{[u°'{x)];x^) 


S'{[u'"{x')]-x'>^) = f^, rx'L'{u'"{x'),u'l,{x')■x'^^) 
= !^d^x det{%){L{u^{x),u^^{x)■x^^) + 6tL} 


(4.3) 


(4.4) 


(4.5) 


Using Jacobi formula for the measure 

f)ST^ 

dV = det(^)d'^x = (1 + ^^)rx, (4.6) 

we get 

6tS{[u‘^{x)];x>^) = J^d'^x{d^6x>^L + 6tL} 

= d’^x{[Ly.]dtU°^ + + ^)dx‘' (4.7) 

+Ma(§^)^tu^-T^jxn}, 

6yS{[u°‘{x)];x>^) = J^d^x{[Lu,.]6yU°‘ + d^{gf^§^6yU°‘)}, (4.8) 

ShSilu^^ix)]] x>^) = J^d'^x{[Ly<.]6hU°‘ + + ^)dx’' 

(4.9) 

= J^d^xd^{L6x>^), 


where [T^q], are the Euler-Lagrange operator and energy-momentum tensor respectively 


[Lu°‘] 

T 

flV 


du°‘ M V duf, ) ’ 


dL 


d{di^u‘^ 


:dyU°' — Lt] 


flU • 


(4.10) 


Thus 6tS = 0, i.e. 6yS = 0 due to Hamilton’s principle and 6hS = 0 due to the invariance of 
re-parameterization of the coordinates that suppose to keep the metric requires to regard 
dyU^" and as independent components and gives rise to the Euler-Lagrange equation 


^Thus the symmetry is the Poincare transformations. In general, the general coordinate transfor¬ 
mations can be considered and all formulae become covariant. 
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(4.11) 


n 

- oA-—) = 0 


du' 


' duZ 


as well as an identity between the Enler-Lagrange operator and conservation property for 
the energy-momentum tensor 


r)T r)T 


(4.12) 


This equation shows that the energy-momentum tensor is conserved if and only if the Enler- 
Lagrange equation is satisfied and the Lagrangian does not depend manifestly on the coor¬ 
dinates. 

The boundary term is as follows 


In + LSx>^)daf^ = 0 . 

If we dehne a current 


: = 


, dL 
du?, 


)6^u°‘ + L6x^, 


(4.13) 


(4.14) 


the equation ([4.13|) leads to an continuity equation for the current if and only if the Enler- 
Lagrange equation is satisfied: 


d 

dxf^ 


= 0 . 


(4.15) 


In fact, this current is just the Noether current with respect to the invariance of the action 
under re-parameterization of the right-angle coordinates on M. 

Remark 4.1: 

Introducing an exterior differential operator dy along the fibre Q M that satisfies 


dy 0, \^dy^ dfi'^ 0, d . dy -I- 


(4.16) 


where dh and d is the nilpotent exterior differential operator on T*X and T*E, E = M x X 
respectively as was as a vertical variational vector held 

^u-=5yU^{x)-^, (4.17) 

then 


5yU°'{x) = i^^dyU‘^ = i^^dE" 


(4.18) 


By means of the vertical variational vector held ([4.17| ) on TQ, 5yS can also be expressed 
as its contraction with 1-form dyS G T*Q 


i^^dyS = 6yS. (4.19) 

We may calculate dyS G T*Q. Since dy commutes with the integral and d"x (see also, for 
example, the functional diherential calculus in [^), it is straightforward to get 
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(4.20) 


d,s= f d^x{[L^c]d,u'^ + ^en}, 

Jn oxf^ 

where 6^ are the Lagrange 1-forms 




(4.21) 


Now by contracting with the vertical variational vector held ( |4.17|) it follows straightforwardly 

s,s. 

Furthermore, due to the nilpotency of it is easy to get 


d 

dySu + — 0 ) 

ox>^ 

where Su is the Euler-Lagrange 1-form dehned by 

£uiu'^ix),u'^{x);x) := 

o;^ are the multi-symplectic 2-forms and in local coordinates: 

02 r n2 r 

UJ^^ ■= = - - -^du^ A du^ + - -—A du^. 


du°‘duu 


du'^dufi 


(4.22) 


(4.23) 


(4.24) 


From the dehnition ([4.23|) , equations ([4.20|) and (|4.22|) , the following theorem |^, [^, [p!4 
holds: 

Theorem 3: For all Lagrangian of a kind of field theories with first order of derivatives 
on the bundle E{X,Q,7i), 

1. The following Euler-Lagrange cohomology is nontrivial: 


Hcft {£u\d£u — 0 }/ {£u\£u — dfi}, 

where [3 = fi{u'^{x),uf{x)]x) is an arbitrary function of {u‘^{x),u‘f{x);x). 

2. The necessary and sufficient condition for conservation of the multi-symplectic 2- 
forms, i.e. 

is the corresponding Euler-Lagrange 1-form being closed. 

Remark 4.2: 

From the dehnition of the Lie derivative it can be seen that the horizontal variations are 
given by the Lie derivative with respect to the variational vector held. 

Let be a vector held on X, exp{Xf) be the how with parameter A, i.e. the one-parameter 
diheomorphism group, induced by ^, / a diherential or a vector on X. The inhnitesimal 
change of / under how is described by its Lie derivative with respect to the vector held f 

LJi^) ■= Jim ^{(/>I/(ea;p(A0a;) - f{x)} = ■^\\=o{filf{x')), x' = exp{Xf)x. (4.26) 
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Here 0^ is the bull-back or the inverse differential for the differential form or vector respec¬ 
tively. 

Taking horizontal variational vector held for the coordinates 


^{x) := 5x^ 


d 

dx^ ’ 


(4.27) 


it follows that the Lie derivative of f{x) with respect to gives rise to the horizontal 
variation of /(x). 

On the other hand, the coordinate variations 6x^ can be expressed by the contraction 
between the variational vector held (4.27) and 1-form dhX^ on T*X, i.e. dx^ on T*E 


i^^dhX^ = 5x^. 


(4.28) 


It is also feasible to express the variation 6hU°‘{x) as contraction of a horizontal variation 
vector held with dyU'^ or dn". To this purpose, ^h,u along the hbre with respect to 
horizontal variations of u°‘{x) 


ih,u ■= 4m"(x) 


du° 


(4.29) 


should be introduced. Combining with the vector held in ([4.27|) , the general horizontal 
variational vector held should be dehned as 


:= + ih,u = 

Its contraction with or du° leads to 

i^^du°‘ = dE^ ■ ih = dhu'^it). 


(4.30) 


(4.31) 


In general, for any functional of m"(x) and m))(x), F(n"(x), n^(x)) : TQ ^ R, its (horizontal) 
variation induced by (ED is 


F(m"(x), m"(x)) F{u°‘{x'),-^U°'{x')) = F(m“(x), m"(x)) -I-dftF(M“(x), m“(x)), 

6hF{u^{x),u'^{x)) = ig^dF(M“(x),M“(x)). 


(4.32) 


Remark 4.3: 

For the total variation, a total variational vector held for m"(x) along the hbre can also 
be introduced 

d d 

^total '■= iv+ih = 

= dx^^— + (d,n“(x) + dx^^n“(x))^, (4.33) 

whose contraction with du" leads to the total variation dtu'^{x) 

ktotaidu"^ = du^ ■ itotai = dtu^{x). (4.34) 
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B. General variation in Hamiltonian formalism 


In order to transfer to the Hamiltonian formalism for classical field theory, we hrst dehne 
a set of “momenta” canonically conjngate to the field variables 


T^a{x) = 


dL 
dii’^ ’ 


and take a Legendre transformation to get the Hamiltonian density 

TTa, VaM“) = 7r„(a;)h“(x) - ft", VaW"), 

where Va = a = — 1. The Hamiltonian then is given by 

H{t) = J^d^-^xH{x), 

with the Legendre transformation 

H{t) = J d'^~^x{7ra{x)u°‘{x) — L{t)}, L{t) = J d^~^xL, 

where S C H is a 3-dimensional simnltaneous space-like hypersnrface in H. 
The action S'([M"(a;)]; ( [4.1|) becomes 

S{[u‘^{x)]-,x^) = [ d'^x{7la{x)u^{x) — H{u°‘,U°‘,'VaU°‘)}. 

Jn 


(4.35) 


(4.36) 


(4.37) 


(4.38) 


(4.39) 


The total variation of the action can be calculated similar to that in the last subsection, but 
here na{x), u°‘{x) and their derivatives should be varied independently. Then we get 


^^^([^"(a;)]; x^) = d’^xld/^Sx'^L -|- StL} 

= d^x{lH^o.]dt7r°^ - 

-h(a^T^^ -F 0r)dx’^ - - T^^Jx’')}, 

where [i7„a], are the canonical operators 

[H,J := u-(x) - £, 

:= *^(x) + - V.(g^). 


(4.40) 


(4.41) 


Similar to the case of Lagrangian formalism, StS = 0, i.e. 6^3 = 0 and 6hS = 0, requires 
to regard 6yU°‘, dv'Xa and Sx^ as independent components and leads to the canonical field 
equations 


= S’ 

an 


'Xa{x) 9^“ + ^«(a(Va«“))’ 


(4.42) 


and an identity between the canonical operators and conservation property for the energy- 
momentum tensor 


[H^.%71'^ - [Hu^]d,u^ + - d,H = 0. (4.43) 

It should be mentioned that all remarks in the last subsection can be made here and the 
theorem 3 can be established as well. 
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V. GENERAL VARIATIONS FOR DISCRETE FIELD THEORY 

We now study the variation problems for the difference discrete held theory with variable 
step-lengths. For simplicity, we consider the cases of 1-1-1 or 2 dimensional hat base manifold, 
i.e. or endowed with suitable signature of the metrics. 

Let be a right-angle lattice on or X^ with nodes X 2 ), n = 1, 2, (i, j) E 

Z X Z and variable step-lengths on 2-directions x^ to be determined by discrete variation 
problems, N be all nodes on L^. For a given node with coordinates Mo := be 

the piece of conhguration space with a set of generic held variables u°‘{x^^’^^) = e Mo 

at the node TM^’d the tangent bundle of with the set of held variables and 

their diherences G the function space on 

Let be the set of nodes neighboring to with index set set of 

nodes related to by the diherences, = [jind{N)\i(i.i) the union of the pieces 

of conhguration space on F(T(Affunction space on T(Af 

Since is a right-angle lattice, it should have only two possibilities for the variable step- 
lengths: either equal step-length variation along two direction simultaneously while along 
each direction the step-lengths are variable, or along one direction the step-length is hxed 
while along the other it is variable. 


A. Variable difference Lagrangian field theory 


The difference Lagrangian for a set of the generic helds m", a = 1, • • •, r, is a functional 
on F(T(M*'®’-^^)) and suppose to be the hrst order of differences of the helds for simplicity 


/^ = 1, 2, 

where as just mentioned, = {xi\x 2 '^), and 


Aiti 


«(*j) — 


u 




(i+1) (i) 

/-Y> ' ' _ /-Y> ' ’ 

X 2 


A 2 M' 


— 


u 




™(i+i) ™(i) 

X 2 ~ X 2 


The discrete action So now reads 

So = y AiiPaA’-Td*’''’, 


(5.1) 


(5.2) 


(5.3) 




where AiXi^ = Xi^^^ 


Xi\ ^ 2 X 2 ^ = X2^^'^ 


AS 

X 2 . 


(5.4) 


Let us consider the coordinates of nodes on the lattice are subject to inhnitesimal defor¬ 
mations that still keep L? as a right-angle lattice 

XiJ) . (AiJ) I xX,j) 

"'"At ^ "'"AtAl M2 ) — “T , 

the corresponding changes in the helds are 

6tu^{xY^’^'> = 6 vU°‘{xY^’^^ -f 5hU^{xY'^’^\ 

6vU°'{xY^’^^ ■= u'°‘{xY^’^^ — u^{xY^’^\ 

6hu'^{xY^’^^ ■= u'^{x'Y^’^^ — u''^{xY^’^^ 

= u°'{x'Y^’^'* — -I- 0(5^) 


(5.5) 
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For the differences of fields, the corresponding changes 


A'^u''^{x'{x' 

W))) 

= A^m° 


i(*j) 

+ (5t A^m 

«(a,)(M), 

5iA^M“(a;)(' 

l,j) = 

5„A^m“ 

(x) 

h.i) 

+ 5hA.M 


(5^A^M“(a;)( 

ij) ■ = 

= A^^i'“( 

xY 

ij) . 



ij) ■ = 

= A>'“( 

x') 

(id) 

- A^m'“( 




= AX 

(x'; 

)(*.i) 

- AXI 



can be calculated to get 


Using the Leibniz law (|3.8|) for differences in each direction, it follows 


5tAiu'^{xY^’^^ = Ai5yU°‘{xY'^’^'^ + ■ AiAi,u°‘{xY'^'^\ 

6tA2U°‘{xY^’^^ = ■ A2Ai,M"(a;)^*’-^A 


The total variation of discrete action ( b-3|) can be calculated 


StSn = E^,J AixP ■ A2xY\AMiLn^^’^'> + A2<54l^(*’^) + 
Using formulae 


Ai 6 x\Ld^''^^ = AYSx\Ld^^-^’^Y - 
A^SxiLo^^’^'^ = A2{64 Ld^^’^-^Y - 


and 




+ 

+ 

+ 

+ 

+ 

+ 






+ acSSSry“ (Ai5a;^(*’^-)) ■ A^m“(*’^')) 

^g^^(A2<5iM“(*’^') - (A2fe^(*’^')) ■ A^n“(*-^')) + 


A ( A ,,a(i- 






we get 


(5.6) 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


(5.11) 
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StSn — SvSd + ShSo 


A ('_9LdA£z 11^X „a(ij) _ ('_£LdA£iA_/^ „«(*j-i) _ (5„ 


aLnC*’-’-!) 


A 


2va(A2w“(*J-i)) 




= J2i,j^iXi^ ■ A2xi^\[L^c(i,j)]6tu‘^^^’^'^ 


■^(1,1'= 




( 6 . 12 ) 


where /i = 1, 2, [L„c(ij)] and are shift operators, discrete Euler-Lagrange operator 

and energy-momentnm tensor respectively 


Ei/h’t) = 


(5.13) 




(i,n : = 


an"hj) 


Ai( 


a(AiM“h-ij))- 


^ a(A2M“(*J“^)) ’ 


(5.14) 




^(AMzz«(*.i)) 




(5.15) 


Regarding and a^-e independent variational bases, ^tS'zj = 0, or d^Sn = 0 

and ShSo = 0, lead to the discrete Euler-Lagrange equation 




— Ao 


'5(A2M“(*4-i)) 


= 0 , 


(5.16) 


a relation between the Euler-Lagrange operator and the (difference) divergence of the discrete 
energy-momentum tensor that may determine the step-lengths 


[L,.(2,2)]A,n“(a;)(*’^) + S^=i,2{A^E-iTz,f+ 




} = 0. 


(5.17) 


It is obvious that all these discrete equation, relation and properties have correct con¬ 
tinuous limits respectively. Furthermore, due to the discrete Lagrangian (|5.1|) depends on 
the differences explicitly, it is possible to introduce the discrete canonical momentum and 
discrete Legendre transformation to transfer to the discrete Hamiltonian formalism as will 
be shown in next subsection. 


Remark 5.1: 

We may introduce exterior differential operators d, dy and dh on T*{M x Xd), T*M and 
T^Xd respectively. They are nilpotent and satisfy 
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d dy dfi^ \^dy^ dfi^ 0. (5.18) 

Especially, dh is due to the difference ow Xd and satisfy Leibniz’s law for ordinary forms 

Remark 5.2: 

Actually, analog to the case with fixed step-lengths , it can be established the dif¬ 

ference version for the Euler-Lagrange cohomology and the necessary and sufficient condition 
for the difference conservation law of the discrete multi-symplectic 2-forms. 


From SySo in (|5.12|) , it is easy to see that we may take dy on Sn to get 

dySo = E A,x\A2xidyLD^^’^\ dyLn^^’^^ = 

(*J) 


(5.19) 


where = 1,2 are the discrete Euler-Lagrange 1-form and multi-symplectic 

potential 1-forms respectively 


y (*,i) ._ rr ,,l^ .MiJ) 

CD •— [-nDya(i,3)\U,yU , 


: = 


dL 


h-ij) 

D 


Then due to the nilpotency of dy, it is straightforward to get 

= 0, := dyeD^^^'^\ 

Therefore, we may get the discrete version for the theorem 3 [^, jl^ : 


orh-i-i) 


(5.20) 

(5.21) 

(5.22) 


Theorem 4-' For all discrete Lagrangian of a kind of discrete field theories with first order 
of differences on the bundle E^Xd^Q^'k) M x Xd, 

1. The following discrete version of the Euler-Lagrange cohomology is nontrivial: 

H]:,pt:={C losed Euler-Lagrange forms}/ [Exact Euler-Lagrange forms}. 

2. The necessary and sufficient condition for conservation of the discrete multi-symplectic 
2-forms, i.e. 


A^iXpCm) = 0, (5.23) 

is the corresponding discrete Euler-Lagrange 1-form being closed. 

Remark 5.3: 

In this paper, L/ is an infinite lattice. It is reasonable to consider a finite lattice. We 
will publish the issues on this topic elsewhere. 


^It is needed some noncommutative differential calculus to completely clarify the properties of 
dh. For the case that Ax^ are fixed, this noncommutative differential calculus can be found in 
|]h|, |R]. For the case of variable step-lengths, similar noncommutative differential calculus can be 
established []T|. 


24 




























B. Variable difference Hamiltonian field theory 


Consider on which there is a right-angle lattice with variable step-lengths in 

each direction, is the base space. 

We hrst dehne a set of the discrete canonical conjngate momenta on the tangent space 
hj)) of the nnion of the pieces of conhguration space on 

which are the set of nodes neighboring to the node 


vr, 


(ti) _ 


(5.24) 


The difference Hamiltonian is introdnced throngh the discrete Legendre transformation 

(5.25) 

The action fnnctional (|0|) now is expressed as 

Sd= Y. (5.26) 

{i,j)£ZxZ 

The total variation of the action can be calcnlated and separated into two parts, i.e. 
the vertical variation SySu and the horizontal variation ShSn 


^tSo — SySf) + ShSf), 

SySo = T,{i,j)ezxz {5yTia^^+^'^^[H^ji+i,j)] - [Hyc{i,j)]SyU°‘^^’^^ 

ShSo = J2{i,j)ezxzAf,x^^^^A{§i^T^J^+^,j)[H^ji+i,j)] - [Hyc(i.j)]5hU°‘^^A 

- E 1 oA ( _L 

^ -‘-Du CX ) 


(5.27) 

(5.28) 

(5.29) 


'^u,ii=1,2{AuEu Tj 




Here 


r ZJ 1_A „,a(i,-i) 

■- AtU - , 


[H^a.{i,3)] :— AiTTo^*’-?’) ~ ^a:( 


a,3) , 


(5.30) 


Regarding 5yU°‘^‘'^\ are independent variational bases, SyS^ = 0 

due to discrete Hamilton’s principle and ShSo = 0 due to discreteized re-parameterization 
invariance on L^, i.e. dtSo = 0, lead to the discrete canonical held equations 


A,7/"h4) = 

aTW+TJT’ 

AiVr^Cr) = 


I A ( \ 

du°‘(X3) “C 


(6.31) 


the canonical form of the relation (|5.17|) that may determine the step-lengths 




} = 0. (5.32) 


It should also be mentioned that all remarks in the last subsection can be made hare and 
the theorem 4 can also be established in this discrete Hamiltonian formalism for held theory. 
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VI. CONCLUDING REMARKS 


In this paper, the variable difference variational approach with variable step-lengths has 
been proposed. It is a generalized version of the difference variational approach with fixed 
step-lengths proposed in |^, [Q. The approach has been applied to both Lagrangian and 
Hamiltonian formalism for discrete mechanics and field theory. Althongh what have been 
dealt with are the systems with first order differences, the key points are available for the 
systems with higher order differences. Obviously, both approaches are different from either 
Lee’s discrete variation with variable time-steps 0,0, m or Veselov’s one with fixed 
time-steps for the discrete classical mechanics [0|, They are also different from the 


discrete variation approach to field theory in that is a generalization of Veselov’s ap¬ 
proach. In our approaches the differences with either variable step-lengths or the fixed ones 
are regarded as discrete derivatives being entire geometric objects. This is more obvious and 
natural from the viewpoint of noncommutative geometry and more analogical to the contin¬ 
uous mechanics and field theory. Therefore, in the continuous limit, the results given here by 
the variable difference variational approach unaffectedly lead to the correct continuous coun¬ 
terparts not only for the equations of motion and symplectic or multisymplectic preserving 
properties, but also for the conservation laws, especially for the energy conservation. 

In view of the structure-preserving criterion for the discrete systems, there are more 
advantages for the variable difference discrete variational approach. Eventually, this has 
been already seen in where in taking the continuous limits for the discrete variation 
problems, which is a generalized version of Lee-Veselov’s variation, combining first discrete 
objects into some difference form is more controllable. 

With variable step-lengths it is, of course, more or less straightforward to generalize 
the symplectic and multisyplectic schemes as ones that are not only symplectic and multi¬ 
symplectic preserving but also discretely energy conserved as has been done for variational 
symplectic energy-momentum integrators in discrete Lagrangian formalism [^, |^, and in 
discrete Hamiltonian formalism [0. But, these discrete formalisms do not transform to each 
other via discrete Legendre transformation. 

The difference variational approach has been applied to the symplectic algorithm and 
multisymplectic one for both Lagrangian and Hamiltonian formalism in |l^. It has been 
shown that the discrete integrants can be combined together in certain manner as certain 
geometric objects in order to construct some numerical schemes with fixed step-lengths as 
variational integrators such as the midpoint scheme in symplectic algorithm and the Preiss- 
man scheme in multisymplectic algorithm. Obviously, the variable difference variational 
approach should be able to apply to the symplectic and multisymplectic algorithms with 
variable step-lengths for both Lagrangian and Hamiltonian formalism. This issue will be 
published in details elsewhere. 

It has been found that the necessary and sufficient conditions for the symplectic 2-form 
preserving in mechanics and the multisymplectic 2-forms preserving in field theories are the 
corresponding Euler-Lagrange 1-form is closed in the relevant Euler-Lagrange cohomology 
il. H- This is also true for the discrete cases |^, [0 as well as the symplectic and 


multisymplectic algorithms |^. For the cases studied in this paper, the Euler-Lagrange 
cohomology should also be true for the various variation problems with variable domains or 
step-lengths. In fact, this matter is already indicated by the boundary terms in the vertical 
parts of these variation problems. We will explain this issue in more detail elsewhere. 
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We have almost completely employed the ordinary description in a coordinate manner in 
this paper not only in order to be more easily understood, especially for non-mathematician, 
but for dealing with both continuous and discrete cases in an analogical manner. It should 
be mentioned however that both the variation problems and the Euler-Lagrange cohomology 
for continuous cases could be dealt with in a coordinate free version in terms of jet bundle 
and variational bicomplex (see, for example, [Q, |P). Although as far as the local issues are 
concerned, the essentials are almost the same. The coordinate free expression should be able 
to contain more general and complicated cases with nontrivial topology. On the other hand, 
however, for the discrete cases the ordinary jet bundle and variational bicomplex approach 
should be generalized to the ones include non-commutative differential calculus in principle. 
We will present the variational bicomplex approach to the issues in this paper elsewhere, 
specially the one with non-commutative differential calculus to the discrete cases 
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